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Abstract 

a: 

Qh! We discuss various aspects of the early petite unification (PUT) of quarks and 

leptons based on the gauge group G PU t = SU (4:) PS ® SU (2) L ® SU (2) R ®SU(2) H . 
This unification takes place at the scale M = 0(1 — 2 TeV) and gives the correct 
value of sin 2 9\y(M^) without the violation of the upper bound on the Kl — > [ie 
rate and the limits on FCNC processes. These properties require the existence of 
three new generations of unconventional quarks and leptons with charges up to 
4/3 (for quarks) and 2 (for leptons) and masses 0(250 GeV) in addition to the 
standard three generations of quarks and leptons. The horizontal group SU(2)h 
connects the standard fermions with the unconventional ones. We work out the 
spontaneous symmetry breaking (SSB) of the gauge group Gptjt down to the SM 
gauge group, generalize the existing one-loop renormalization group (RG) analysis 
to the two-loop level including the contributions of Higgs scalars and Yukawa 
couplings, and demonstrate that the presence of three new generations of heavy 
unconventional quarks and leptons with masses 0(250 GeV) is consistent with 
astrophysical constraints. The NLO and Higgs contributions to the RG analysis 
are significant while the Yukawa contributions can be neglected. 



1 Introduction 

The idea of Grand Unification (GUT) based on simple groups like SU(5) [DI2] or SO(fO) 
[HUE], characterized by a single gauge coupling (/gut, is a very attractive scenario of the 
physics beyond the Standard Model (SM). In GUT models quarks and leptons are gen- 
erally members of the same representation under the given gauge group, and this results 
in transitions that violate quark and lepton quantum numbers. As such transitions are 
very suppressed in nature, the GUT scale must be very large, typically C(10 16 GeV), in 
order to be consistent with the experimental data, in particular with the lower bound 
on the proton life-time. 

A less ambitious program is the Petite Unification 0E] that aims at unifying quarks 
and leptons at some energy scale M, not much greater than the electroweak scale, with 
the gauge group Gs <8> Gw that is characterized by two independent couplings gs and 
gw- It is further assumed that Gs and Gw are either simple or pseudosimple (a direct 
product of simple groups with identical couplings). An attractive choice for the strong 
group Gs is SU(4) a la Pati-Salam [7j with the lepton number playing the role of the 
fourth colour. It turns out jH] that with this choice of Gs only very few weak groups 
Gw can have low unification scale being simultaneously consistent with the measured 
value of sin 2 6 l i^(M|), the upper bound on the rare decay Kl — > lie and the data on 
flavour changing neutral current (FCNC) processes. Basically only two weak gauge 
groups SU(2) 3 and SU(3) 2 can be made consistent with the experimental data. 

The general properties of the S77(4)ps ® SU(2) 3 and S77(4)ps <8> SU(3) 2 unifications 
have been discussed in jSJ. In these models the values of sin 2 9w at the unification scale 
M turn out to be 1/3 and 3/8, but a very fast renormalization group evolution allows 
to obtain correct sin 2 6 w{M^) with M « 1 TeV and M m 3 TeV, respectively. 

Concentrating on the unification based on the gauge group 

Gput = ST/(4)ps ® SU(2) L ® SU(2) R ® SU(2) H , (1.1) 

let us recall three most interesting properties of this model: 

• In addition to the standard three generations of quarks and leptons, new three gen- 
erations of unconventional quarks and leptons with charges up to 4/3 (for quarks) 



and 2 (for leptons) and masses 0(250 GeV) are automatically present. The hor- 
izontal group SU(2)h connects the standard fermions with the unconventional 
ones. 

• The placement of the ordinary quarks and leptons in the fundamental representa- 
tion of SU(4)ps is such that there are no tree-level transitions between ordinary 
quarks and leptons mediated by the S77(4)ps gauge bosons. This prevents rare 
decays such as Kl — > \xe from acquiring large rates, even when the masses of these 
gauge bosons are in the few TeV's range. 

• There are new contributions to flavour changing neutral current processes (FCNC) 
involving standard quarks and leptons that are mediated by the horizontal SU(2)h 
weak gauge bosons and the new unconventional quarks and leptons. However, they 
appear first at the one-loop level and can be made consistent with the existing 
experimental bounds. 

It should be emphasized that the existence of new heavy fermions that are placed in 
the fundamental representation of S77(4)ps together with ordinary fermions is essential 
for having early unification of quarks and leptons without the problems with the rare 
decays like K L — > fie. The original petite unification group SU(4)p$ eg) SU(2) 4 jH] having 
only ordinary fermions in the fundamental representation of S77(4)ps and consequently 
Kl — > fie proceeding at the tree level is ruled out as the early unification model unless 
further new physics, such as large extra dimensions, is invoked jH]. In comparison with 
the usual left-right symmetric models based on the S , L r (4) PS ® SU{2) L ® SU{2) R group 
Ell, that is relevant for the grand 50(10) unification, the present model has an 
additional SU(2)h factor and the fermion representations that differ from the latter 
case. These new ingredients allow a low unification scale that is not possible in the 
SU(A) PS ® SU[2) L ® SU[2) R model. 

In the present paper we would like to extend the analysis of the unification in (jl.lj) 
presented in || by 

• working out the spontaneous symmetry breaking (SSB) of the gauge group Gput 
down to the SM gauge group, 

• generalizing the one-loop renormalization group analysis of [6 J to the two-loop level 
and including the contributions of Higgs scalars and Yukawa couplings, 



• demonstrating that the presence of three new generations of heavy unconventional 
quarks and leptons with masses (9(250 GeV) is consistent with astrophysical con- 
straints. 

A short discussion of the rare decay Kl — > /ue and of FCNC processes was already 
presented in jH] and will be elaborated on elsewhere. 

Our paper is organized as follows. In Section 2 we recall the main ingredients of 
the model in question, presenting in particular the fermion representations. In Section 

3 we present the Higgs system that accomplishes the desired SSB of Gput down to the 
SM group and work out the formulae for gauge boson and fermion masses. In Section 

4 we set up two-loop renormalization group equations for the evolution of the gauge 
couplings and determine the petite unification scale M using as the inputs sin 2 W (M§), 
a s (M§) and a(M§). We present a very simple formula for M as a function of these three 
inputs that to a very high accuracy reproduces our numerical analysis. In Section 5 we 
address the fate of the new heavy fermions in the context of astrophysical constraints. 
Our conclusions and a brief outlook are given in Section 6. 

2 The Model 

2.1 Preliminaries 

In this section we will describe the main ingredients of the model based on the group 
Gput m (jl-lj) - After presenting the pattern of the SSB and of the conditions for the 
coupling constants, we will describe in turn the gauge boson sector and the fermion 
sector. The Higgs system responsible for the SSB will be presented in Section 3. 

2.2 General Structure 

The pattern of the SSB breaking is assumed to be 

Gput ^ G x ^ G 2 ¥$ SU(3) C ® U(1) EM , (2.1) 

where 

G ± = SU(3) c (g 3 ) ® U(l) s (~9s) ® SU(2) L (g 2L ) ® SU(2) R (g 2R ) ® SU(2) H {g 2H ) , (2.2) 

and 

G 2 = SU(3) c (g 3 ) ® SU{2) L {g 2 ) ® U(l) Y ( 9l ) (2.3) 



is the SM group. In order to streamline the notation we will denote the usual hypercharge 
coupling g' by g% and (?qcd by #3, reserving the index "5"' for the U(l)s group. Moreover, 
if not specified, the coupling g 2 will always stand for g 2 L- 

Thus at some scale M, the strong S77(4)ps group is broken down to the product of 
the gauge symmetry group of QCD and the strong U(l)$ group, corresponding to the 
unbroken diagonal generator T 15 of the S77(4)ps group that does not belong to SU(3) C . 
The explicit expression for T 15 is given in (J2.2UJ) . 

In the next step, at scale M < M, the subgroup U(l) s ® SU(2) R ® SU(2) H of the 
G\ group is broken down to the weak hypercharge U(l)y group. The generator Ty of 
U(l)y is given by 

Ty = C R T 3R + ChT 3 h + CsTi5 (2.4) 

with T 3R and T 3 # being the diagonal generators of SU(2) R and SU(2)h, respectively. 
Consequently, the electric charge generator is given by 

Q = T 3L + T Y = Q W + C S T 15 (2.5) 

where Qw is the "weak" charge corresponding to the group Gw 

The coefficients Cj in (J2.4J) describe the embedding of the weak hypercharge U(l)y 
group into G\. The fact that the weak hypercharge U(l)y group merges into both G$ 
and Gw at M, allows us to put quarks and leptons into identical representations of the 
weak group Gw and consequently make the quarks and leptons to be indistinguishable 
when the strong interactions are turned off. In the model in question, the coefficients Cj 
are given by 

C W = C 2 R + C 2 H = 2, Cl = | (2.6) 

These two values play an important role in the RG analysis presented in Section 4. 
In particular, C w appears in the crucial group-theoretical prefactor 

*** = T*!r = \< (27) 

which, for the present discussion, is simply 1/3. The reader is urged to consult |5 a and 
jH] for prefactors corresponding to other choices of Gw- 

Renormalization group effects in the range Mz < yU < M can decrease the mixing 
angle down to the experimental value sin 2 ^(Mf ) w 0.23, provided the unification scale 
and the representations of the matter fields (fermions and scalars) are properly chosen. 
An explicit one-loop relation between sin 2 9w(M z ) and sin 2 6 W is given by [2J E] 

sin 2 ^(M 2 ) = sin 2 C[l-^^7^-8vr-a(M 2 )(irin^- + K'ln^)], (2.8) 



where 

a (Ml).^l, a 3 (Mi),iM), (2.9) 

and the coefficients K and K' are given in Section 4. 

It has been demonstrated in [H] that with the values of Cfy and C| in (j2.fi |) and 
the fermion representations specified below, the resulting values of K and K', allow to 
obtain the correct value of sin 2 #vi/(M§) provided M ~ M ~ 0(1 TeV). At the two-loop 
level that we discuss in Section 4, the RG analysis is much more involved and we will 
proceed differently. We will use the experimental value of sin 2 0^(M|) as an input to the 
evolution of the gauge couplings and will determine the unification scale M, in analogy 
to GUT analyses, by studying the "matching" conditions for the relevant couplings. 
While these conditions will be spelled out systematically in Section 4, the basic relations 
behind them are [U] 

+ r t^t? . ( 2 - 10 ) 



e 2 (Mf) [g 2 (M 2 z W MM 2 )] 2 ' 
92l(M 2 ) = g 2R (M 2 ) = 92h(M 2 ) = 9w(M 2 ) , (2.11) 



W +7^4^ (2-12) 



[ft(M 2 )] 2 [gw(M 2 )] 2 [g s (M 2 )} 



with g w = g 2 and 



g 3 (M 2 ) = ~g s (M 2 ) = g s (M 2 ) . (2.13) 

Finally, we will use the standard MS definition for sin 2 9\y, namely 

sin 2 9 W {M 2 Z ) = 4^11- (2-14) 



9!l(M, 



2.3 Gauge Bosons 

2.3.1 Gluons and Leptoquarks 

The SU(3) C content of the adjoint representation 15 of SU(4:)ps is as follows: 

15 = 8 + 3 + 3 + 1, (2.15) 

and the corresponding gauge fields A 1 ^ are represented by 



/ 



{^(1 = 1, ...,15)} 



.4 



ttS 



(z = l,...,; 



Gp'\ 



Gp 2 



Gp 3 



KG- 1 G-i "-■< 



(2.16) 



G^ ' A^s J 



Here the octet A^ s stands for the gluons, 



r<±A 






Gr 



—ft { A tis =F tA fis) 
77= { A fis =F lA/is) 
-ft { A tis T iA.^) 



(2.17) 
(2.18) 
(2.19) 



A]^ s is the neutral gauge boson which corresponds to the generator T15 of 



and A^ s 

SU(4)s and equivalently to the generator of U(l)s- Explicitly 



T, 



15 



2^ 



/l 



V 



\ 



-3/ 



(2.20) 



The leptoquark gauge bosons Gp l carry electric charges ±| and connect quarks to 
leptons. They are responsible for the rare transitions, like Kl — > yue, the phenomenology 
of which has been briefly presented in jHj and will be presented in detail elsewhere. Under 
the breaking of SU(4)s down to SU(3) C x U(l)s the leptoquarks G^ 1 gain masses of 
order M whereas the gluons and the gauge boson A^s remain massless. 

2.3.2 Electroweak Gauge Bosons 

The model has nine massive electroweak gauge bosons in addition to the massless photon. 
These include (i) four charged gauge bosons Wj^, W^ and two neutral gauge bosons 
Z\i Z 2 all with masses of order M, and (ii) the standard W and Z° gauge bosons with 
the conventional masses that are very precisely measured. Below we summarize the 
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structure of the neutral gauge boson sector. The derivation of these formulae and the 
discussion of the gauge boson masses is postponed to Section 3. 

It should be remarked that the field B^ of the SM is expressed in terms of the fields 
A^s and (W^j , which couple to the diagonal generators of U(1)s®SU(2)r(S)SU(2) h , 
respectively, as follows: 

B, = cos 9 S A, S + ^ (W* R + W^) , (2.21) 

where the mixing angle 9s is defined by 

tan0<j = ^-— . (2.22) 

9w 2 



Furthermore, we have 



~ cos q 

Z x » = sin 9 S A, S - — yj- (Wl R + Wl H ) , (2.23) 

Z 2 , = -L (wl H - W^ R ) . (2.24) 

The masses of Z 1/Lt and Z 2 ^ are related through the relation 

M Zl = -^f-, (2.25) 

COS Oq 



which will be derived in Section 3. Equation ()2.25|) is the analog of the SM relation 
Mw = Mz cos 9w M must be larger than 800 GeV in order for the model to be consistent 
with the experimental data. Finally, it follows from (J2.12JI that the hypercharge U(l)y 
coupling constant g\ is defined in terms of gs and gw by 

VSgsQw 9w sin e s 

(6g s + 8 9w) v2 

with all couplings evaluated at \x — M . This equation is analogous to the well-known 
relation e = gsm6 w . 

2.4 Fermions 

The fermions in each generation can be divided into two groups (with the electric charges 
shown in parentheses): 



a) Ordinary Fermions 

with ^ and ^ transforming as (2, 1) and (1, 2) under SU{2) L <g> SU(2) R , respectively. 

b) New Heavy Fermions 



( ^(4/3) ^ 

L,R \" ay "' / L,R 



Ql,R= A' > W= r J > ( 2 - 28 ) 



V^(i/3)y 




with unconventional charges and masses of 0(Mp)- The left-handed and right-handed 
fields transform as (2, 1) and (1, 2) under SU(2) L <g> SU(2) R , respectively. Thus their 
transformation properties under SU(2)l ® SU(2) R are the same as of ordinary fermions 
but due to different electric charges, they cannot be considered as new generations of 
ordinary quarks and leptons. The three generations of these new heavy fermions con- 
stitute for themselves a new set of sequential generations, that are connected with the 
ordinary generations through the interactions mediated by the SU{2)h gauge bosons 
and the leptoquarks Gv' 1 . 

In jHj also the third group of vector-like fermions 

,,„ _ ( U'-"(5/S) \ _ /(~S»(-l/2)\ 

with unconventional charges and masses of O(M) has been considered. Here Q' L R and 
L' L R transform as doublets under SU{2) L and are singlets under SU{2) R ® SU(2) H , 
whereas Q" L R and L" L R transform as doublets under SU{2) R and are singlets under 
SU (2) l® SU (2) h ■ These additional fermions relevant in principle for the renormalization 
group analysis for scales \i> M were introduced in [U] to keep the SU(2) L , SU(2) R and 
SU{2) H couplings equal under renormalization group evolution at least at the one loop 
level. Meanwhile we have realized that this equality is broken already at the one loop 
level by Higgs contributions and consequently there is really no reason to introduce them 
at all. Moreover, as discussed in Section 5, due to very strange charges the existence of 
these fermions is problematic for cosmology. Therefore we will exclude them from our 
analysis. 

Until now we have given only the representations under SU(2)i and SU(2) R . In order 
to avoid dangerous tree level transitions involving ordinary fermions that are mediated by 



the S77(4)ps gauge bosons and FCNC transitions mediated by the "horizontal" SU(2)h 
gauge bosons, we proceed as follows: 

• With respect to S , f/(4) PS , we put ordinary quarks ip q L R together with the heavy 
leptons Llr into the fundamental representations of this group and similarly for 
ordinary leptons i/) l L R and heavy quarks Ql,r- 

• With respect to SU(2) H , we put ordinary quarks ip q L R together with heavy quarks 
Ql,r i n the doublets of this group and similarly for ordinary leptons ip l L R and 
heavy leptons Ll,r- 

Explicitly then, the fermions are placed in the representations under Gput as follows: 

/ (rf c (l/3),f>(4/3)) (lJ-l),v(0))\ , s 

^ L = (4,2,1,2 ) L = v y i h \ i a i^i h \ )> \ 2>30 

- V K(-2/3), D(l/3)) U-2),l(-l))) L 

( (rf c (l/3), UU/3)) (U-l),u(0))\ , s 

* fl =(4,l,2,2) fl = V Kl ' y n \ K ! V !! \ (2.31) 

\(u%-2/3),D(l/3)) U-2),l(-l)) J r 

where in order to put the ordinary quarks into representations of Gput we used their 
charge conjugated fields, suppressing a minus sign in front of the u c field. Each entry in 
the 2x2 matrices in (J2.30J) and (J2.31)) represents an SU(2)h doublet, whereas the columns 
represent SU(2)l^r doublets. With respect to 5*?7(4)ps we do not show explicitly the 
colour indices but <i c (l/3) is placed in the quartet together with l u (— 1), £7(4/3) together 
with i/(0) and analogously for fermions in second rows in ()2.30|) and (12.311) . 

3 The Spontaneous Symmetry Breaking 

3.1 Choices of Higgs Scalars 

In this section, we will discuss the Higgs sector which is needed to spontaneously break 
SU(4)p S <g> SU(2) L <g> SU(2) R <g> SU{2) H down to the SM group. In particular, we will 
focus on those scalars that can make an important contribution to the RG evolution of 
the gauge couplings of the SM. 

The Higgs scalars that are needed are the following: 

• The Higgs field that breaks SU(4) PS can be written as 

$ s = (15, 1,1,1). (3.1) 



The vacuum expectation value (VEV) of this Higgs field breaks <S77(4)ps down to 
SU(3) C <g> U(l) s at a scale M. 

Below M, one has effectively SU(3) C <g> 5C/(2) L <g> 5C/(2) H <g> 5C/(2) H ® £/(l)<?. One 
would like to find a scalar field that can spontaneously break SU(2)l <g> SU(2) R ® 
SU(2) H <S> U(1)s down to SU(2) L ® U(1) Y at some scale M > M z . Such a scalar 
should transform non-trivially under SU(2) R ® SU{2) H <g> £/(l)s, in particular it 
should have a non vanishing U(l)s quantum number. Furthermore, we require the 
U{\)y gauge boson, B^, to be a linear combination of W^ R , W^ H , and A m5 , which 
implies that there should be mixing between these three gauge bosons. On first 
look, it appears that there are several possibilities. 

One might consider, for example, the following Higgs Fields: 

A" = (4,1,3,1) A" =(4,1,1,3), (3.2) 

where a — 1, ..4 denotes the SU(4)ps index. 

For symmetry reasons one might want to include the following Higgs field: A^ = 
(4,3, 1, 1). However, its VEV would break SU(2)l and there exist severe experi- 
mental bounds on the contribution of any Higgs triplet to the p parameter (already 
at tree level). These bounds imply that its VEV would have to be less than a few 
percent of the SM VEV. For convenience, we will assume that this VEV is identi- 
cally zero. 

Finally there are Higgs fields which would break SU (2) l®U (1)y down to U(1)em- 
Furthermore, these Higgs fields should also couple to fermions as in the SM so as 
to give these fermions a mass. The fermions in our model transform as: (4, 2, 1, 2)^ 
and (4, 1, 2, 2) R . A mass term would transform as the bilinear 



(4, 2, 1, 2) L <g> (4, 1, 2, 2) R = (1 + 15, 2, 2, 1 + 3). (3.3) 

Therefore, in principle, we could have the following Higgs fields: (15,2,2,1), 
(15,2,2,3), (1,2,2,1), (1,2,2,3). Any one of these Higgs fields can be a suit- 
able candidate for the symmetry breaking. The question one might ask is whether 
or not they are all necessary. To study this question, let us first concentrate on 
the S'?7(4)p S -singlet Higgs fields. We consider first 

$ = (1,2,2,1). (3.4) 

10 



As we shall see explicitely below, the VEV of $ gives equal masses to all fermions 
(quarks and leptons, including both conventionally and unconventionally charged 
fermions). This fact alone necessitates additional Higgs fields. Let us then look at 

$^=(1,2,2,3). (3.5) 

The VEV of this S'C/(4)pg-singlet, SU (2) jj-triplet Higgs field would split the mass 
scales of the conventional fermions (both quarks and leptons) from the unconven- 
tional ones, as shown below. This, however, does not split the masses of the quarks 
and the leptons. For this to happen, we again need additional Higgs fields. 

To split quark and lepton masses (for both conventional and unconventional fermions), 
it appears sufficient to just use a Higgs field: 

<\P = (15, 2, 2,1), (3.6) 

where (5 — 1, .., 15. As we shall see below, in the coupling of (ft 13 to fermions, it is 
convenient to write 

0, = /^, (3.7) 

where A^/2 are the generators of S77(4) PS . Notice that, under the QCD subgroup 
SU(3) C , a 15 splits into 8 + 3 + 3 + 1. Therefore it should be the SU(3) C singlet 
part which develops a VEV. A VEV of the form (<f) s ) = (0 15 )Ais/2 would split 
the masses of the quarks from those of the leptons. The origin of further splitting 
in the quark sector, in particular between the top quark and the other standard 
quarks will be discussed in Section 3.3. 

Finally, for symmetry reasons, one might also have (1, 2, 1, 2), (1, 1, 2, 2) and (15, 2, 1, 2) 
Higgs representations. However, these scalars do not couple to fermions. The electric and 
hypercharge structures of these Higgs fields are identical to those of the representations 
(1, 2, 2, 1), the color-singlet part of (15, 1, 2, 2), and (15, 2, 2, 1) respectively and we will 
concentrate on the latter. 

3.2 Electric Charges of Selected Higgs Fields and their Cou- 
plings to SM Gauge Bosons 

Since we are primarily interested in this paper in the contributions of the scalar sector 
to the RG evolution of the SM gauge couplings up to M, we shall not discuss the case 
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of $5 = (15, 1, 1, 1). We shall instead concentrate on A Rt n, $, &h, 4>s- All Higgs fields 
in this paper are complex. 

As we have seen in jHUHj and above, the U(l)s charge matrix for SU(4)p,s fundamental 
representations is 

fl/3 \ 

1/3 

1/3 

\ 0-1/ 

We now use (13. 8J) to find the electric charge assignments and the U{1)y quantum numbers 
for the Higgs fields listed above. 



Ys — C S T 15 



(3.8) 



For both A R and Ah, one has 



Q W = (1,0,-V 



(3.9) 



since they are both triplets under their respective Gw group. 

With Q = Qw + Ys, one obtains the following electric charge assignments for both 

A r ,h- 

Q A = {(4/3, 1/3, -2/3), (0, -1, -2)} , (3.10) 

where the first entries inside the parentheses are for the color triplets and the 
second entries are for the color singlets. 

One can now explicitely write A^ and A# . One has: 



U/3),i a (1/3),* A (-2/3),iN / A o 



Ar, h - {(A RH ' , A RH '\ A RH '*), (A RH , A RH , A RH )} 



(3.11) 



where i = 1, 2, 3 is the color index. 

Since A R and A H are SU{2) L singlets, the U(l)y quantum numbers Y are identical 
to Qa, namely 

Y = Q A . (3.12) 

In the discussion of symmetry breaking given below, it is convenient to express A R 
and Ah as 2 x 2 matrices, namely 



Arw = A R)H -Tr,h 



^R^/y/2 -AW2 



(3.13) 
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for the color singlet part, and 



Ag^/2 Ag»-/V2 



A- RH = A' R _ H .T RM - [ A M/3,, i/V 5 _ A <w /2 I • < 3 - 14 > 

for the co/or triplet part. The superscript 5 in (J3.13J) is a notation for "color- 
singlet". For symmetry breaking, only the neutral scalar in (J3.13J) can develop a 

VEV. 

For $ = (1,2,2,1), it is clear that one now has two SU(2)l doublets. Therefore 
each doublet interacts with the SU(2) L gauge bosons in a standard way. As we 
have already shown in jB], the charge assignment for a representation (1, 2, 2, 1) is 
simply 

Q = Qw = f _" 3 i M , (3.15) 

where the first row (column) refers to T 3 ^ = 1/2 (T 3R = —1/2) and the second row 
(column) refers to T 3 l = —1/2 (T 3 # = 1/2). 

One can explicitely write $ in terms of a 2 x 2 matrix as 

$ = 1 2 . 3.16 

The U(l)y quantum number is in this case simply 

Y = T 3R = ±^, (3.17) 

where Y = —1/2 for the SU(2)l doublet with charge (0, —1) and Y = +1/2 for 
the SU{2)l doublet with charge (1,0). 

The charge assignment for <3># = (1, 2, 2, 3) is just slightly more complicated. It is 
simply a direct sum of Qw (I3.15|) and T 3 # = (1, 0, —1) (for a triplet), namely 

Q = Qw®T 3H . (3.18) 

From (|3.18|) . one obtains the following 2x2 charge matrices referring to (2,2) 
under SU{2) L ® SU{2) R and there are three of those (with the U(l)y quantum 
numbers listed next to them): 



H )iy=(i,+5), (3.19) 
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(3.20) 



(3.21) 

In these equations, the U(l)y hypercharges Y refer to the first and second column 
respectively. Also equations (|3.19|) . (|3.2(J|) . (|3.21|) refer to six SU{2) L doublets 
which couple to the corresponding gauge bosons in a standard way. The couplings 
to the U{1)y gauge boson are given in terms of the hypercharges listed above. 

Finally, one can write explicitely $# as 

3*, + =(!f + ' 1 I! + ' 2 ), (3-22) 




*h,o = ( "^ J"'' ) , (3.23) 

$„_ = f^I' 1 ^' 2 ] , (3.21) 

\®H-,1 ®H-,2 

The subscripts +, 0, — in the equations above refer to the three scalars with T%h = 
(1,0,-1) respectively. 

The Higgs field, whose VEV could split the masses of the quarks from those of 
the leptons, could be (ft 13 = (15, 2, 2, 1) as mentioned above. To find the charge 
structure of this Higgs field, it is useful to first split it into SU(3) C representations: 

1)8: 0*(i = l,..,8), 

2) 3,3: ± ' 1 = ^(0 9 T # 10 ), ± ' 2 = ^(0 n T *0 12 ), ± ' 3 = ^(0 13 T # 14 ), 

3) 1: 15 . 

Notice that each one of the 0's written above is a (2, 2) under SU(2)l ® SU(2)r. 

The £/(l)s quantum numbers, Yg, of 0^ can easily be found. They are exactly the 
same as the electric charges of the 5 , [/(4) P s gauge bosons for the latter are singlets 
under Gw- Therefore: 

1) Y s = for 8, and 1, 
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2) Y s = ±| for 3,3. 

Since Q = Qw ® Ys, one obtains: 




(3.25) 



4/3 7/3 \ 5 11 

1/3 4/3/ ' V 6' ' l ' 

■4/3 -1/3 \ 11 5 

-7/3 -4/3 i 6' 6' 



Explicitely, we can write 0^ as 



;(8) 



,« 



01 


4»>^ 


01 


4 8, '° , 


rfw 


4*+N 


4 1 ''" 


4*» , 


,(3),4/3 


*f-" 3 


,(3),l/3 


4 3),4/3 



,(3) -4/3 ,(3) -1/3 



(3.28) 



(3.29) 



(3.30) 



^ (3) -^k-7/3^3),-4/3j • ( 3 - 31 ) 

From the above listing of the SU(2)l doublets and their hypercharge quantum num- 
bers, one can proceed to include their contributions to the evolution of the SM gauge 
couplings at one loop. At two loops, in order to properly include the scalar contributions, 
one has to work out the Yukawa couplings between some of the Higgs fields listed above 
and the fermions. 

Finally, as we mentioned above, in principle there could exist, for symmetry reasons, 
the following Higgs fields: (1, 2, 1, 2), (1, 1, 2, 2), and (15, 2, 1, 2). The electric charge and 
hypercharge assignments for these Higgs fields are identical to those for $, the 1 of </r , 
and <f)@ respectively. 
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3.3 Yukawa Couplings 

In this section, we will not make any serious attempt to construct a model for fermion 
masses, but rather we are more interested in a rough value for the Yukawa couplings as 
deduced from the overall mass scales of the conventional and unconventional quarks and 
leptons. This exercise serves as an estimate of the contributions of the Yukawa couplings 
to the two-loop beta functions. For convenience, we shall make two assumptions con- 
cerning the masses of the conventional and unconventional fermions. We believe that our 
estimates of the unification mass scales will not be much affected by details of fermion 
mass models. These assumptions are as follows. 

1) What we have in mind for the discussion that follows is some kind of democratic- 
type mass matrices for the ordinary quarks ^T] . This implies that the masses for the Up 
and Down sectors given below will be taken to be universal mass scales which appear in 
front of 3 x 3 matrices whose elements have magnitudes of order unity. For this reason, 
models of this type are usually referred to as Universal Strength for Yukawa couplings 
(USY) because these couplings are common to all families for each Up and Down sector. 
The hierarchy in masses comes from the diagonalization of this type of matrices. The 
above ansatz could be realized in a number of ways. One approach is to go to more 
than three spatial dimensions as has been done in |T2] where it was found that the quark 
mass matrices were of the almost-pure phase type (a generalized version of democratic 
matrices). In J3], this type of mass matrices was found to fit the mass spectrum and 
the CKM matrix rather well. 

2) For the unconventional fermions, we have to keep in mind that they are as yet 
unobserved. As a result, their masses are constrained in various ways which, in general, 
depend on their lifetimes and decay modes. If, for the sake of making some estimates as 
to their contributions to the RG evolution, we assume that their masses are all equal to 
or larger than 250 GeV, then we are faced with a very different mass pattern for these 
unconventional fermions. The 3x3 matrices mentioned above for the ordinary quarks 
would have to be very different for the unconventional ones. In fact, most likely they 
would have to be of a form as to yield eigenvalues of the same order so as to generate 
masses which would differ from each other by at most a factor of two. 

With the above two assumptions, the discussion which follows deals uniquely with 
universal mass scales. 

The fermions of our model are for each generation as follows: 

^ L = (4, 2, 1, 2) L , V R = (4, 1, 2, 2) R . (3.32) 
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As we have discussed at length in |>| and in Section 2, these representations contain 
conventionally and unconventionally charged quarks and leptons. As alluded to above, 
one needs to split the masses of the conventional fermions from the unconventional ones 
since the latter have not been observed experimentally This can be achieved by the use 
of both $ and $ H . 

Then, the quark and leptonic SU(2)h doublets are 

(m^R, Ql,r), (Ll,r^ 1 l , r ). (3.33) 

We emphasize again that the way the fermions are ordered here implies that there is 
no tree-level transition between normal quarks and normal leptons due to the SU(4)p$ 
gauge bosons which link only conventional to unconventional fermions. 

With (J3.33J) in mind one can now try to write down the various Yukawa couplings. 
First, it is useful to note the following identity: 



(zr#)(^r) = (^L^r- (3.34) 

Since $# is an SU(2)h triplet, it is convenient to write 

®h-T h = „ „ , r- . . , (3.35) 

where $^" = ($^ T i& H )/V2. 

Using (|3.34|) . one obtains for the Yukawa coupling to $: 

g^L^R + h.c. = g*{(il>i L $>!)* + Ql$Qr 

+ (quarks — > leptons)} + h.c. , (3.36) 

where g$ is taken to be positive. 

Similarly, the Yukawa coupling to $# can be written as 

- g* H V L Q H -f H * R + h.c. = g^ H {Q L ^Q R - (ip q L -^VrY + (quarks 

— > leptons) + (charge — changing 
interactions) + h.c.} , (3.37) 

where g^, H is taken to be positive. 
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When $ and $^ develop a VEV, (|3.36|) and 13.37)) give contributions to fermion 
masses. First let us remind ourselves that both $ and $^ are 2x2 matrices with 
respect to SU(2)i ® SU{2)r. Therefore we have a different VEV for the Up and Down 
quark sector (as well as for its leptonic counterpart): 



(3.38) 
(3.39) 



o ($) d 
(*2r>« 

We list below various contributions to fermion masses. 
I) Masses coming from $ and $#. 

• The Up quark sector: 

m u (&,& H )=g*(&)u-g* H (®H)u, ( 3 - 4 °) 

• The Down quark sector: 

m D {$, $ H ) = g*($) d - g$ H ($ 3 H ) d , (3.42) 

mo(Q 1 $H)=g*($)d + 9* B (®H)d- (3-43) 

For reasons which are given below, we will (by a suitable definition of the phase) 
define the lepton masses with a negative sign ( which cannot be detected in any 
case). We have the following masses. 

• The Up lepton sector: 

- m Zu ($, $ H ) = g$(®) u - gt> H (® 3 H ) u , (3.44) 

- m N ($, $ H ) = 0*($>u + g$ H (®h)u ■ (3.45) 
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• The Down lepton sector: 

-m Ld ($,$ H ) = g*(*U- 9* H (*H)d, (3-46) 

-m L (Q,$ H ) =g${®) d + g<!> H {® 3 H)d- (3.47) 

II) Quark-lepton mass splitting from 0^ = (15, 2, 2, 1). 

As we have mentioned above, the component of 0^ = (15, 2, 2, 1) which can develop a 
VEV is the color singlet <jP~\ most specifically cpi and (fi 2 in ()3.29|) . Since (p^ = 15 , 
it follows that 

(<^7» = (0 15 )T 15 , (3.48) 

where diagonal matrix Ti 5 is given in ()2.20|) . Let us denote the Yukawa coupling of the 
fermions to </r by g^. 

The total contribution of the fermion masses is now given by the following expressions. 
Notice that we only give the expression for the neutral lepton sector as Dirac masses. 
No attempt will be made in this paper concerning possible nature of the neutrino mass. 
This will be dealt with in a subsequent paper. 

• The Up quark and lepton sector: 

mu = g*(®)u - g* H {®V)u + gMi ], °) , ( 3 -49) 

m = g 9 (Q) v + g$ H (& H ) u + gMi )fi ) , (3-50) 

- m Lu = g*(<S>) u - g* H (®H)u ~ 3^(^ 1) '°> , (3.51) 

- m N = g^) u + g* B {& H ) n - 3g4<j>? fi ) . (3.52) 

• The Down quark and lepton sector: 

m D = gi($) d - g* H ($ 3 H )d + 94>(<f>2 )fi ) , (3-53) 
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(1),0\ 



m D = g*($) d + g* H (& H ) d + gM h ) , (3-54) 



- m ld = g^) d - g* B (& H ) d - 3^(#'°> , (3.55) 

- m L = g*($) d + g* H ($ 3 H ) d - 3^(0^'°) . (3.56) 

From (J3.49J) . (J3.50J1 . ()3.51|h one obtains the following products of the Yukawa couplings 
with the VEVs: 



1 
1< 



0*(*) u = -{mu + 2rrifj - m L ) , (3.57) 



9*h($hU = \( 2m u ~ 2m u) , (3-58) 

gM )fi ) = \(rn Lu + mu ), (3.59) 

9^}d=j(mD + 2mf)-m Ld ), (3.60) 

9* B {&H)d = \(2m 6 -2m D ), (3.61) 

9M )fi ) = \(^L d + ^D). (3.62) 

From the above equations, one can estimate the size of various Yukawa couplings by 
making educated guesses on the masses of various unconventional fermions as well as the 
values of various VEVs. 



3.4 The Neutral Gauge Bosons from the Breaking SU(2)r (g) 
SU(2)h®U(1)s^>U(1)y 

The main goal of this subsection is to find the expression for the U(l)y gauge boson B^ 
in terms of the neutral gauge bosons of SU(2) R eg) SU(2) H <g> U(l)s- However, we will 
also present similar expressions for the other two massive neutral gauge bosons which 
accompany the massless B^ at this stage of symmetry breaking. We will postpone the 
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discussion of mixings between charged gauge bosons to Section 5 where we discuss the 
fate of the unconventional fermions. 

Before discussing how the neutral gauge bosons obtain their masses through sponta- 
neous symmetry breaking, let us explain how (|2.21jl was derived. From (|2.4|) for Ty, one 
can readily write down an expression for the U(l)y gauge field as 



B u = 9w = C S A„ S + 9S = V C iW W ui . (3.63) 

V n 2 (-12 _|_ % 2 r 2 ^ n 2 O 2 Z n 2 O 2 ^— ' " M K ' 

9w L/ s "+■ 9s^w V SW^s "t y5 b w i^su{2) L 



Putting into (J3.63J) the explicit values C| = 8/3, C^ = 2, one obtains for the coefficients 
in the notation of (j2.21|) explicitely 

a 2 9w . a \3gs ( Q RA s 

C ° S 0s = /o~2 , a 2 = ' Sm ^ S = /o~2 I A 2 Z ■ ( 3 - 64 ) 

Furthermore, from (|2.12j) . it is straightforward to derive ()2.26|) . 

Let us now see how one obtains (j2.21|) and the two massive neutral gauge bosons by 
explicit coupling to the Higgs fields. 

In this discussion, it is convenient to write the SU{2)r,h gauge bosons in terms of a 
2x2 matrix as follows: 

W? - = -( W " R ^ ^ W ^A (365) 

(R ' H) '" 2 2 \V2W^ -W^ H) J- 

As mentioned above, the Higgs fields that can accomplish this breaking should carry 

the quantum numbers of SU(2)r (g) SU(2)n <S> U(l)s- They are Ar = (4, 1,3, 1) and 

Ah = (4, 1, 1,3). The components of these Higgs fields which can acquire a VEV are 

the following color-singlet and neutral scalars: A RH . From (J3.13J) . one finds: 

(a|h) ^o<a W/v ^W^. (366) 

The next step is to calculate the masses of the neutral gauge bosons. This is shown 
in Appendix A. Here we just quote the results. 

The mass matrix squared for the three neutral gauge bosons is found in Appendix A 
to be 

1 9 2 w 5l -y/S/Sgwgs^R N 

g 2 w 6 H -VW*9wgs6 H ■ (3.67) 

\ y/3/89wgspR -VmgwgsSjj (3/8)~g 2 s (5 R + 8%) j 
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It is straightforward to diagonalize this matrix. The expressions for the general results 
are rather long and are given in Appendix A. Here we just present a special case which 
is quite reasonable in its own right, namely 



ir = o H 



5. (3.68) 



With (J3.68|) and Appendix A, we obtain the following mass eigenstates and eigenvalues. 

B, = cos 9 S A^s + ^ {Wl R + Wfa), M B = 0, (3.69) 

where cos 9 s and sin 9 s are defined in (|3.64|) . 
2) Z llt : 



Z llx = sin 9 S A, s - ^ (WJr + W*b), M Zl = \ ^3^ + 4^5. (3.70) 

3) Z^: 

Z 2ll = -L(W; H -W* R ), M Z2 =g w 8. (3.71) 

From (J3.70J) and (J3.71J1 . one obtains also the interesting relationship in (J2.25J1 . 

Apart from a slight difference in expressions and notation, the above presentation is 
very similar to the one given in j3] for the group SU(4)ps <S> SU(2) A . 



4 Renormalization Group Analysis 

4.1 Preliminaries 

In this section we will generalize the renormalization group (RG) analysis of the gauge 
couplings and of sin 2 9w presented in jU] to the two-loop level. Moreover we will include 
the contributions of the Higgs scalars to the relevant /3-functions that were neglected 
there except for the standard Higgs doublet. At the two- loop level the scalars contribute 
to the running of the gauge couplings both via gauge as well as Yukawa coupling. The 
prime goal of this analysis is to check whether the inclusion of these new effects does 
not spoil the early unification of quark and leptons analyzed at one-loop level in j^j. In 
calculating the relevant j3 functions we used the general formulae of [T4*j . 
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We will consider, as in j^j, two scenarios. One with M = M and the other with 
M > M. In both scenarios we will set the masses of new fermions in (J2.28J) as well as 
those of the Higgs fields (|3.4jl . (J3.5J) and (|3.fij) to be equal to a single scale Mp with 

M F = M 0i$H> ^ = (250 ± 50) GeV, (4.1) 

while we will assume that the scalars in (J3.2J) have masses very close to M so that their 
contributions to the gauge couplings evolution for renormalization scales fi < M can be 
to first approximation neglected. On the other hand, these contributions must be taken 
into account for scales /i in the range M < \x < M. 

In both scenarios, that is M = M and M > M, we will use as the experimental 
inputs the values of 

ai{M z ) = — — a 2 {M z ) = . 2 . 2 , , a 3 {M z ) 4.2 

cos z u w {M z ) sm W (M Z ) 

with the MS values [T3] 

l/a(M|) = 127.934(27), a 3 (M|) = 0.1172(20), (4.3) 

sin 2 9 w (M 2 z )\ cxp = 0.23113(15) . (4.4) 

The couplings a,i will then be evolved by means of the RG equations presented below 
to the scale Mp at which the contributions of the new heavy fermions (J2.28J) and of the 
Higgs fields (|3.4|) . ()3.5|) and (|3.6JI to the (3 functions will be switched on. Again, in order 
to streamline the notation we will denote the usual a 1 by a\ and ckqcd by «3, reserving 
the index "S"' for the U(l)s group. 

In the scenario with M = M, the unification scale is simply found from the petite 
unification condition (J2.12J) 

- W +7-^77, (4-5) 



MM 2 )] MM 2 )] MM 2 )] 

with Cw and C| given in (|2.6|) . In what follows it is useful to denote the l.h.s of (J4.5J) 
by Fp and the r.h.s. by Fr, that is 

i r 2 r 2 

and F R = -^- + -%-. (4.6) 



[«i(M 2 )] [« 2 (M 2 )] [« 3 (M 2 )] 

In obtaining (|4.5jl we have used 6>s(M 2 ) = a^{M 2 ) that is valid for M = M. The 
condition ()4.5|) corresponds to the equality of the properly normalized gauge couplings in 
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Figure 1: Fl and Fr in (|4.6|) as functions of the renormalization scale /i for central input 
variables. The unification scale for M = M is uniquely determined by the crossing point. 
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Figure 2: a 3 and as as functions of \x > M for central input variables. The unification 
scale for M > M is uniquely determined by the crossing point. 
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the standard GUTs as SU(5) and 5*0 (10). Graphically the scale M is found as depicted 
in Figure ^ where we plot Fl and Fr as functions of the scale \x. The crossing point of 
these two evolutions determines uniquely the petite unification scale M = M. 

In the scenario with M > M, the evolution of a\, a 2 and 0:3 for scales jjl < M is as 
in the scenario with M = M, but now the condition (J4.5J) is replaced by 

W + 7^77, (4-7) 



MM 2 )] [a 2 (M 2 )] [a s (M 2 )] ' 

in accordance with (|2.28jl . with as corresponding to the U(l)s group. Formula (|4.7jl 
allows to determine the value of as{M 2 ) at a given scale M that should be larger than 
800 GeV in order to be consistent with the lower bound on the right-handed gauge boson 
masses. 

Above M the gauge symmetry group is SU(3) C ® U(V)s ® SU{2) 3 and the evolution 
of the gauge couplings must also include the contributions of the Higgs scalars (|3.2|) . 
The unification scale M is then simply found from the condition 

a 3 (M 2 ) = a s (M 2 ). (4.8) 

This is illustrated graphically in Figure El As at the two-loop level the evolutions of the 
couplings 0:3 and as are affected by the three SU(2) gauge couplings we also need their 
values at M. This is simply found from 

a 2R (M 2 ) = a 2H (M 2 ) = a 2L (M 2 ) (4.9) 

with a 2 ^(M 2 ) = a 2 (M 2 ) calculated as in the scenario with M = M. 

In what follows we will first give the RG equations relevant for the case M = M. 
Subsequently we will study the scenario with M > M. 

4.2 Renormalization Group Equations (M = M) 

4.2.1 The Range M z < fi < M F 

The relevant RG equations for the evolution of the couplings a\, a 2 and a% are given as 
follows 



dai a 2 

dfi 2n 



5=1,2,3 



(4.10) 



with i = 1,2,3 and A* = g 2 /47r, with g t being the Yukawa coupling of the top quark. 
We neglect the contributions from Yukawa couplings of lighter quarks and from the 
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standard leptons. Moreover as the evolution of the Yukawa couplings is rather involved 
at the two-loop level, we will keep them at constant values corresponding to m t (m t ) in 
the case of the top quark with the same procedure for new heavy quarks and leptons 
discussed below. This turns out to be a very good approximation in our case, where the 
evolution of couplings takes place over a rather short range of scales, but of course such 
a procedure could not be justified in the case of GUTs. 

The coefficients (/?o)i, (Pi)ij an d {Pi)Jt are weu known but one has to remember that 
the U(l)y coupling a,\ used here is differently normalized than the a% coupling in the 
SU(5) model. With three generations of quarks and leptons we have then 



(A))i = -y, 



(A)) 2 = y, 



(A)) 



7. 



(4.11) 



0x 



Z-95/9 -3 -44/3 \ 

-1 -11/3 -12 
V —11/6 -9/2 26 



and 



ml 



17 

y 



(/3i)l 



3 
2' 



Wl 



(4.12) 



(4.13) 



4.2.2 The Range M F < ^ < M 

Above the scale Mp the contributions of new fermions in (|2.28jl and the scalars in (J3.4J) . 
(13. 5 j) and (J3.6J) have to be taken into account. This modifies the coefficients in (J4.11J) 
and (J4.12J) as follows 



(A 



136 



on 



209 



(A))s 



2 
3 



19 

y 



(A) 



16 

y 



(4.14) 



/ -2230/9 - 17777/27 -42 - 209 -280/3 - 3200/9 ^ 
$i= -14-209/3 -158/3-247/3 -24-128 , (4.15) 

y -35/3-400/9 -9-48 -50-592/3 j 

where second terms in (I4.14|) and (J4.15J1 represent scalar contributions. In addition 
Yukawa couplings of new quarks and leptons have to be taken into account. With 
respect to the democratic fermion mass model presented in Section I3~31 the Yukawa term 
in the RG equations (|4.10|) is generalized to 



A y /i 



doti 
djji 



2tt 



£w 



Y *k 

lAfc 4^' 



(4.16) 
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where the sum runs over 13 Yukawa couplings corresponding to seven heavy quarks and 
six heavy leptons. The coefficients {(3i)J t are as in (J4.13J1 . The coefficients {(3i)J k with 
k 7^ t relevant for the contributions of new fermions are found to be 

• with respect to SU(3) C : 



(ft)a3 - 2' 




(&)& = \> 




m Y lQ = 


89 


(00* = 


29 


mXi = 


13 

T 


(A)Jk = 


25 

Y 



(A)*q = 2, Q = U,D,S,C,B,T (4.17) 

with respect to SU(2)l' 

Q = U,D,S,C,B,T 

L = l u ,ld,ls,lc,h,h (4-18) 

with respect to U(1)y'- 

sq . . . 

Q = U,C,T 

Q = D,S,B 

L = l d J s J b . (4.19) 

It is clear from these formulae that the Yukawa couplings of (U,C,T), (D,S,B), 
(l u , l c , l t ) and (Zd, l s , /(,), generally differ from each other. However, in order to get a rough 
estimate of Yukawa contributions let us set all Yukawa couplings of the new fermions to 
be equal to the top Yukawa coupling g t . Then effectively only the X t term is present as 
in (J4.1UJ1 but the coefficients (Pi)J t are modified as follows 

71 3 97 

ml = -^ ml = y , ml = ^- c^o) 

4.2.3 Analytical Formula for M 

We will present the numerical analysis of this scenario in subsection 4.4. On the other 
hand it is possible to derive an approximate analytical formula for M as function of 
the input parameters 1)4.1)1 . (J4.3)) and (J4.4J) . To this end we define effective "one loop" 
coefficients 

/3f = Wi + E (A)«S + (/ ^Y' (4 ' 21) 

j=l,2,3 
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with olj and At frozen at some intermediate scale M z < \x < Mp. Analogous procedure 
is used for the range Mp < fx < M = M. 
We then find 



M = M F 



M z lKa/Kt "' 



exp(P), (4.22) 



where 



and 



r2 o(M)_ sin 2 9 W (M 2 ZJ . 

°S„. /^2\ "' „:_2/)0 J ^""J 



8na(M 2 )K tot V «s(M|) sin 2 # 



if 



^3 = ^ (/?f CI + /3f ^ - /3f ) (4.24) 

with the RG coefficients given in (J4.11|) - (j4.12|) . Formula (J4.24J) is also valid for K tot but 
this time the coefficients (|4.14j) - (j4.2U|) relevant for the range Mp < {i < M should be 
used. The formula (J4.22J) can be directly obtained from ()2.8|) by setting M = M and 
making the replacement: 

M M F M 

Kin ► K 3 In — + K tot In —. (4.25) 

M z M z M F K ' 

4.3 Renormalization Group Equations (M > M ) 

4.3.1 Preliminaries 

The evolution of couplings from /i = M z to \x = M proceeds as in the previous scenario 
and consequently the formulae given in the subsections 4.1 and 4.2 allow the determina- 
tion of the couplings a 3 (M 2 ), a s {M 2 ), a 2 L{M 2 ), a 2 R(M 2 ) and a 2 H{M 2 ), that constitute 
the starting point for the subsequent evolution from M to M. We already mentioned 
that the Higgs fields break the equality of the three SU(2) couplings and consequently 
five different couplings have to be considered, although the splitting of the three SU(2) 
couplings is insignificant as one can see in Figure El 

4.3.2 The Range M < fi < M 

The relevant RG equations for the evolution of the couplings a 3 (/i 2 ), as([i 2 ), «2l(/^ 2 ), 
«2_r(/^ 2 ) and a 2 .f/(/i 2 ) are given as follows {i = 3, S, 2L, 2R, 2H) 

** = (**\ +A Y (^), (4.26) 

where the first term on the r.h.s describes the evolution of the gauge couplings with all 
Yukawa couplings set to zero and the second term takes into account the presence of 
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Figure 3: Running of 05, a^, o^h for M > M, with all input variables taken at their 
central values. 



these couplings. We then have 

da. 
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dn 



G 



a 2 



(A)« + £09 



Dij 



4:71 



(4.27) 



with j running over the five couplings and 



*(*©--£»«£. 



(4.28) 



as in flUHjj) but this time i=3, S, 2L, 2R, 2H. 

With three generations of ordinary and heavy quarks and leptons the coefficients in 
(J4.27JI are given as follows (second terms below stand for scalar contributions): 
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and with the ordering i = 3, S, 2L, 2R, 2H 



A "< A UH A V 



(4.29) 



(4.30) 
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where 

/-50-1 -9 -9 -18 \ 

-8 -7 -9 -9 -18 

(&) = -24 -3 -158/3 -18 (4.31) 

-24 -3 -158/3 -18 

v -48-6 -18 -18 -452/3 j 

are the contributions arising from fermions and gauge bosons and 

/ -658/3 -67/6 -48 -60 -12 \ 
-268/3 -277/6 -48 -60 -12 
(A) = -128 -16 -247/3 -57 -24 (4.32) 

-160 -20 -57 -157 -24 
^ -32 -4 -24 -24 -448/3 j 

are contributions coming from the scalars. With the separation of the scalar contributions 
the breaking of the symmetry between the SU(2) couplings becomes obvious: ol2h-,olil 
and a 2 R evolve slightly differently. For the coefficients (A)ir we find 

• with respect to SU(3) C : 



(Pi)I Q 



2, Q — all quarks 



(4.33) 



with respect to SU(2)l- 
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(4.34) 



(4.35) 



with respect to S77(2)#: 



(A)L„ = 3, 



(0 

with respect to U(l)s- 
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(4.36) 



(4.37) 
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with all remaining coefficients being zero. 

Again as in the case of the range Mp < \x < M we can get a rough estimate of the 
effects of the Yukawa couplings by setting them equal to each other. In this case 

,2 



where Ay is universal and 



Ay [fj, 



Ws 



da.i 

dfi 



a i (a \Y^Y 



(4.38) 
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(4.39) 
(4.40) 



4.3.3 Analytic Formula for M 

Proceeding as in subsection 4.2 one can derive an approximate analytic formula for M 
as a function of M, M F and the input parameters (J4.1J1 . (|4.3|l and (|4.4|) . We find 



M = M 



~M Z ~ 

_M F _ 


K 3 /K' 


~M F 
.M . 



Ktot/K' 



exp(P') 



where 



P> = p 



K tot 
K> 



with K-i, Ktr<i and P defined in subsection 4.2 and 



K' 



CI 



167T S 



(Pf'Pf) 



(4.41) 



(4.42) 



(4.43) 



where /3| ff and /3| ff are evaluated as in (J4.21J) but with various coefficients relevant for 
the range M < \x < M. 

4.4 Numerical Analysis 

We have solved the RG equations listed above numerically to find M = M in the first 
scenario and M in the scenario with M > M and M > 800 GeV. Let us first neglect 
the Yukawa contributions. The dependence of M on the input parameters is in the first 
case as follows: 

• As sin 2 Qw(Mz) is very precisely known, the variation of its value within the full 
range in (|4.4|) introduces only a shift in the ballpark of 30 GeV. M decreases with 
increasing sin 2 6V(^f )• 
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• The uncertainty of M due to Mp and a^ (M§) is significant as shown in tabled 
where we have set a(M|) and sin 2 #yi/(M§) at their central values and varied Mp 
and az(Mz) in the ranges given in (|4.1jl and (J4.3JI . M increases with increasing 
M F and a 3 (Mf). 

• The effect of the NLO contributions is significant. They increase the scale M = M 
by roughly 160 GeV 

In the case M > M our findings are as follows: 

• As expected, the value of M increases with decreasing M. In table Q the last two 
columns show the LO and NLO values of M for M = 800 GeV. The pattern of the 
dependence on Mp and 0:3 (M§) is similar to the one in the case M = M. 

• The NLO corrections in this case are even more important and amount to the 
increase of M by 370 ± 40 GeV for M = 800 GeV with a smaller shift for M > 
800 GeV. 

Most importantly 

• The inspection of the renormalization group coefficients of Sections 4.2 and 4.3 
shows that the large impact of NLO corrections originates to a large extent in the 
scalar contributions. 

• Moreover, as discussed below, the latter contributions increase M roughly by 
250 GeV and 450 GeV for M = M and M > M = 800 GeV, respectively. 

In summary, varying all the input parameters within one standard deviation we find 
at the NLO level 

M = (1253 ± 132) GeV, (M = M) (4.44) 

and 

M= (1713 ±302) GeV, (M > M = 800 GeV). (4.45) 

4.5 Anatomy of Renormalization Group Effects 

In what follows we would like to present the results of a number of exercises that we 
made in the context of our numerical analysis. 

First the inclusion of Yukawa couplings changes the values in (I4.44J1 and (J4.45J) re- 
spectively to 

M = (1273 ± 132) GeV, (M = M) (4.46) 
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Table 1: Values for the unification scale M for various values of 03 (M§) and Mp, l/a(M|) = 
127.934 and sin 2 6 W {M 2 Z ) = 0.23113). In the last two columns we set M = 800 GeV. The 
RGEs have been set up including gauge, fermionic and scalar contributions. 



«s(Mi) 


M F 


M = M, LO 


M = M, NLO 


M> M, LO 


M> M, NLO 




200 


961.4 


1120.9 


1083.0 


1410.4 


0.1152 


250 


1036.7 


1180.8 


1226.2 


1538.6 




300 


1102.6 


1233.7 


1357.2 


1655.6 




200 


1012.1 


1190.4 


1178.5 


1562.3 


0.1172 


250 


1091.3 


1253.0 


1334.5 


1701.6 




300 


1160.7 


1308.2 


1477.0 


1828.9 




200 


1063.6 


1262.3 


1279.0 


1726.4 


0.1192 


250 


1146.9 


1327.3 


1448.2 


1877.3 




300 


1219.7 


1385.0 


1602.9 


2015.2 



and 

M = (1742 ± 300) GeV, (M > M = 800 GeV) . (4.47) 

The relatively small impact of Yukawa couplings on our results justifies our rough treat- 
ment of these contributions. 

Next, removing the scalar contributions altogether we find instead of (J4.44J1 and (J4.45J) 



and 



M = (994 ± 120) GeV, (M = M) 



M = (1254 ± 300) GeV, (M > M = 800 GeV), 



(4.48) 



(4.49) 



implying that these contributions are very important. 

Finally, we have investigated the impact of the increase of Mp to M = M. This 
is equivalent to the removal of the heavy fermion contributions to the renormalization 
group equations below the unification scale. We find then for the maximal unification 
scale obtained with largest a 3 (M|) and smallest sin 2 6 l vy(M|) 



M max = M = M F = 4595 (2442) GeV 



(4.50) 



with and (without) scalar contributions (at the scale 250 GeV), respectively. Thus 
even without heavy fermions below /i = M we obtain a relatively low unification scale. 
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That is, the group theoretic factors Cs and Cw in our master formula guarantee early 
unification even if only the SM fields have masses below the unification scale. On the 
other hand the presence of new heavy fermions with masses O (250 GeV) is necessary in 
order to keep M below 2 TeV. 

5 The Fate of the Unconventional Fermions 

5.1 Relevant Interactions 

As we have shown in [Hj and in Section 2.4, the construction of the model necessitates 
the introduction of new heavy quarks and leptons with unconventional charges: ±4/3 for 
the quarks and ±2 for the leptons. These fermions behave in exactly the same manner 
as the ordinary fermions and are linked to the latter by SU(2)h gauge interactions 
9w^ l,bWii,h-Th 7 m ^ l,Ri and by the Yukawa interactions g<s> H ^f l$ 'h -Th^ r. The question 
arises as to the fate of the lightest of these new fermions. 

One should note that the vector-like quarks and leptons as written down in [6 a and 
(J2.29J) suffer from problems with cosmological constraints in the following way. Since they 
carry electric charges such as 5/6 or 3/2, the lightest one would be absolutely stable since 
it cannot decay into known particles. There are severe constraints on objects such as 
stable fractionally charged leptons [T7] and, unless the numbers are severely depleted by 
exotic mechanisms such as the one proposed in [T7], they are ruled out observationally. 
We will therefore omit them altogether. 

Let us first discuss the fate of the unconventional SU(2)h partners of standard quarks 
and leptons by looking only at the gauge interactions. For the lightest new quark, if 
the gauge bosons do not mix, one would face the dreadful conclusion that it would be 
absolutely stable. We will discuss below various constraints which rule out that situation. 
Fortunately, we will see that there are Higgs fields whose VEVs mix the gauge bosons of 
the different groups and it is this mixing which renders the lightest new quark unstable, 
i.e. it can now decay into conventional fermions which are assumed to be lighter. A 
similar argument applies to the case of the new leptons. 

In this section, we focus only on the charge-changing interactions and hence on the 
charged gauge bosons of SU(2)i (g> SU(2)r <g> SU{2)h- We have already discussed above 
the breaking of SU(2)r ® SU(2)h and, in particular, the neutral gauge boson sector. 
In that discussion, we have used the Higgs fields A# = (4, 1, 3, 1) and A# = (4, 1, 1, 3). 
However, as we have shown in Sections 13.11 and 13.21 other Higgs fields come into play, 
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namely $ = (1, 2, 2, 1), $# = (1, 2, 2, 3) as well as ft 1 ' which is the color singlet part of 
0^ = (15,2,2, 1). The VEVs of these scalars will evidently mix various gauge bosons. 
We have also mentioned earlier that, for symmetry reasons, one would also like to have 
A L = (4,3,1,1), <f>i = (1,2,1,2) and $ 2 = (1,1,2,2). A detailed treatment of this 
problem is given in Appendix A. Here we will give a simplified version for illustration. 
In particular, we will make the same assumption as in ()3.68|) . namely 5r = 5h = 5. 
Furthermore, we will assume that 5 > ($ 2 ) ^ (3*), (&h), (^i); (<P }• F° r simplicity, 
we will also assume that (A/,) = 0. Below, we will denote a generic electroweak scale by 
v and a generic large scale by S. 

Since this section focuses primarily on the question of whether or not the lightest 
of the unconventional fermions are stable, we will present a streamlined version of the 
discussion which will accurately summarize the points that we wish to make. Some of 
the details can be found in Appendix A. 

The lightest of the unconventional fermions can only decay into either a conventional 
fermion plus the SM W boson or into three conventional fermions if it has the appropriate 
mass. This implies that the current J 1 ^ should interact with the SM W for the former 
case or it mixes with either J^ or J^ for the latter case. This is what we will show below. 
We use the notation gw for the three gauge couplings which are defined to be equal to 
each other at M (see (12. lip . In the estimate of the decay rates given below, we shall, 
however, take the value of g 2 or equivalently of the Fermi constant Gp for simplicity. 

With the above remarks taken into account, we now list the eigenvalues and mass 
eigenstates for the charged gauge bosons. (The exact expressions with various factors 
included are given in Appendix A.) Our notations are as follows. We use Wl,Wi,W 2 
for the mass eigenstates and Wl, Wr, Wh for the gauge eigenstates. We have 



M Wl = 0{g w v), M Wl = 0{g w 6), M m = 0{g w 5). 
The relationship between the gauge and mass eigenstates is given by 



(5.1) 






o 



T 



(wt\ 



I i 

0(v 2 /6 2 ) 
k O(v 2 /6 2 ) 



0(v 2 /S 2 ) 0{v 2 /5 2 )\ 



i 

1 

V2 



1 
~V2 ) 



(wt\ 



(5.2) 



Let us denote generically the charge-changing currents associated with SU(2)l, SU(2)r, 
and SU(2)h by J£, J^, and J^ respectively. The charged current interactions involving 
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these currents can now be written as (with ± omitted for simplicity) 

^-interaction = 9W\y Li "Ri ^ HI * " fJ.R 

\W,h) 

In terms of the gauge boson mass eigenstates, one can now write (15. 3|) as 

fw tlL \ 

^interaction 9W\" Li Ri HI 



(5.3) 






(5.4) 



where 



(5.5) 



i^Li Jri J H) ~ i^Li J Ri ^ H) X ^ ' 

and where the rotation matrix T is given in ()5.2|) . Explicitely, we now list separately 
the following interactions. 



Interaction with W^l- 



£l = 9w Jl w ^l 

= g w [r L + 0(v 2 /5 2 )(J^ + J^}W, L 



(5.6) 



Interaction with W^i'. 



£i = 9w JrWui 



9w[0(v 2 /5 2 )J"l + ^(J"r + J'hW,i 



Interaction with W^'- 



£>2 = 9W Jh W »2 

= g w [0(v 2 /6 2 )J L l + -L(J£ - J H )]W, 2 



A few remarks are in order here. 



(5.7) 



{hi 



From (J5.6J) . we observe that the contribution of a V + A current to known weak 
interactions is suppressed in the interaction Lagrangian by a factor 0(v 2 /6 2 ). Ex- 
perimental constraints [15] from searches for right-handed currents in normal weak 
interactions give 0(v 2 /5 2 ) < 10~ 3 . This can easily be satisfied within our model 
through the choices of the various VEVs. 
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• Also from (|5.6jl . one can see that 3^, which changes a conventional fermion into 
an unconventional one and vice versa, couples with W^l with a factor 0(v 2 /S 2 ). 
This has some important implications concerning the decay of the lightest of the 
unconventional fermions. 

1) If the lightest unconventional fermion which appears in 3^ has a mass greater 
than the sum of the accompanying conventional fermion mass and the W-boson 
mass, it can have the decay mode F — > f + W (W being the standard W^l) via 
the interaction 

g w O(v 2 /5 2 )3£W flL . (5.9) 

Even though 0(v 2 /5 2 ) < 1CT 3 , the decay rate can be substantial because the 
unconventional fermion decays into a real Wl- 

2) If the mass of the lightest unconventional fermion is less than the sum of the 
VT-mass and the mass of the accompanying conventional fermion, the decay can 
occur through the interaction: 

2g 2 w O(v 2 /S 2 ) 31 - X 2 3, H . (5.10) 

1 ~ m w L 

• If the mass of the unconventional fermion is between the sum of the accompanying 
conventional fermion mass and the W-boson mass and either that of W± or W 2) 
the decay would be into a real Wl as in (1) above. 

• If the mass of the lightest unconventional fermion is less than the sum of the W- 
mass and the mass of the accompanying conventional fermion then, from (J5.7J1 and 
(J5.8J) . one obtains the following additional contributions to the decay: 

1 X 

u q 2 - m 2 wl q 2 - m 2 W2 

In computing the matrix elements for various decays, one has to express the fermions 
which appear in 3£, 3^ and 3^ in terms of the mass eigenstates. This will result in the 
appearance of a number of mixing angles which are different from the CKM elements. 
In the absence of a plausible model of fermion masses (especially for the unconventional 
ones) , the best one can do is to make estimates of the decay rates based on reasonable 
assumptions about the magnitudes of these unknown mixing angles. 
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5.2 Decay Modes 

We will present here some estimates of possible decay modes of the lightest unconven- 
tional quark and lepton. The main purpose, as we have mentioned above, is to see how 
fast or how slow these fermions decay. Since there are strong constraints on "stable" 
fermions from cosmology and from collider experiments, the unconventional fermions 
should be sufficiently heavy and should decay fast enough to evade these bounds. This 
is what we will show below. A comprehensive study of all possible decays is beyond 
the scope of this paper and it will be included in a future publication. Here we wish to 
merely present some illustrative examples. 

To be more specific, let us assume that the lightest unconventional quark is £7(4/3) 
and the lightest unconventional lepton is /«(— 1), where the numbers inside the paren- 
theses represent the charges of these particles. (One can easily change this scenario to, 
e.g., _D(l/3) and ld(—2), or any other combination.) Since those lightest unconventional 
fermions cannot decay into other unconventional fermions, the only option left is for 
them to decay into conventional fermions. 

In what follows, we will use the following notations for the conventional fermions: 
ifj d = d,s,b and ipl = z/ e , z/ M , u T , for the quarks and the neutral leptons respectively. 
Notice that SU(2) H interactions connect £7(4/3) to ip d , and /„(— 1) to ip l v . 

In using J 1 ^, we will express the fields which appear in that current in terms of their 
mass eigenstates. Consequently, the results presented below will contain mixing angles 
of the type V 0ipy 

• £7(4/3) decay: 

We now consider the following possibility: M(£7(4/3)) > m{ip d ) + M^ . This case 
is more or less obvious since we require the unconventional fermions to be heavy, 
i.e. around 250 GeV or so. 

In this case, the dominant decay mode would be the semi-weak process: 

£7(4/3) - ft + W+ L . (5.12) 

Since, by assumption, m dsh /M'~ J <C 1, one can immediately find the following decay 
width: 

G F M 3 - M\ 2 M\ 
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where Gf is the Fermi constant. For Mjy r* 2bQGeV ) one obtains 

r^6|^^| 2 |(9( W 2 A 2 )| 2 GeV, (5.14) 

where V^« is the matrix element of V = U^Ufj. The matrices C/d and ^ are 
those that diagonalize the mass matrices of the conventional Down-quark sector 
and the unconventional Up-quark sector, respectively. 

The mean lifetime is found to be 

r » 1.1 x io- 25 |v^r 2 i^v^-v (5.i5) 

For illustration, let us put 0(v 2 /5 2 ) ~ 1CT 3 as we have mentioned above. The 
mean lifetime is r^ r; 1.1 x 1 19 1 V^ ,» | -2 s. One can see that U can decay very 
fast unless the mixing Vjjm is abnormally small. 

For a particle which decays that fast, there is no cosmological constraint. It can 
be searched for at future facilities such as the Large Hadron Collider. This type of 
search was discussed at length in [To] . 

• l u (—l) decay: 

The computation for the decay rate here is very similar to that presented above. 
We will assume that the mass of l u (—l) is comparable to that of U. The main 
decay mode is then 

l u (-l)^^ + W- L , (5.16) 

where tp l u = u e , v^ v T . 

Once again we will assume that M~ L rj 250 GeV. The decay width is then found 
to be 

r L *G\V Iur f\0(v 2 /5 2 )\ 2 GeV. (5.17) 

Its mean lifetime could be very short, i.e. r^ w 1.1 x 10 — 19 | V^ ,, |~ 2 s, if \Vj u ^i | 2 is 
not abnormally small. Again there is no cosmological constraint. A discussion of 
the search for leptons of this type can be found in jl6j. 

As we have mentioned above, although we have chosen U and l u to illustrate how 
an unconventional fermion can decay entirely into conventional particles, one can choose 
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other unconventional fermions to be the lightest ones and study their decays. This will 
be presented elsewhere. The main point in this section was to show that, because of 
mixing among the various gauge bosons, the lightest unconventional quark or lepton is 
fairly unstable, and, for the range of masses that we consider, decays mainly into a real 
W and a conventional fermion. 

6 Summary 

In this paper we have extended the discussion of the early unification of quarks and 
leptons based on the gauge group SU(4) PS ®SU(2) L ®SU(2) R ®SU(2) H jU]. In particular 

• we have presented the Higgs system which accomplishes the spontaneous symmetry 
breaking (SSB) of the gauge group Gput down to the SU(3) C (g) £7(1)qed group 
with the acceptable spectrum of gauge boson, fermion and Higgs masses. 

• we have shown that the inclusion of NLO effects and of Higgs scalars into the 
renormalization group analysis increases the unification scale M, relatively to the 
estimates in [6j, by roughly 250 GeV and 450 GeV for M = M and M > M = 
800 GeV, respectively. This allows still for a unification of quarks and leptons at 
scales 0(1 — 2 TeV). Specifically in the two scenario considered we find 

M = (1253 ± 132) GeV, (M = M) (6.1) 

and 

M= (1713 ±302) GeV, (M > M > 800 GeV). (6.2) 

• we have shown that the presence of three new generations of heavy unconventional 
quarks and leptons with masses 0(250 GeV) is consistent with constraints coming 
from cosmology. 

A detailed discussion of the rare decay Kl — > [ie and of FCNC processes will be presented 
elsewhere. 
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A Appendix A 

In this Appendix we present the mixing of gauge bosons following the symmetry breaking 
of U(l) s ® SU(2) L <g> SU(2) R <g> SU(2) H -► U(l) Y <g> SU(2) L -► C/(1) QED . 



A.l Neutral gauge boson mixing 

The mixing of neutral gauge bosons W 3R , W 3 h, A$ is given by the VEV of A R (4, 1, 3, 1 
and A#(4, 1, 1, 3) (considering only color singlet part) 

'0S R /y/2 




(Ah) 



<Ah> 



'0 6 H /V2* 




with the corresponding covariant derivatives 



r 



D M A fi = d^A R - V3/8 ffsA^(-l)Afl - ig w [w; R —, A R 



T" 



D,A H = 0„A H - V3/8 ^(-1)^ - Ww[W; H -, S H 



and 



" 2 2 ^W- -W 



(A.l) 
(A.2) 

(A.3) 
(A.4) 

(A.5) 



/i 3 



From Tr ( D fl A R D^A R J and Tr ( D^A h D^Ah ) follows the squared mass matrix of 
neutral gauge bosons 



/ 



9w s % 







g 2 w 5 2 H -^/W9w~9s5 2 h 

\-VW8gw~g s 5 R -^/mgw~g s 5 2 H f<?i(<% + <%) / 

This matrix can be diagonalized by an orthogonal matrix O: 

-VWgw~g s s R \ /o o o \ 

o <?^<% -y/z/SgwgspH 

\-V3/%9w~gs5 2 R -Vmgw~gs5 2 H l~g 2 s (K + K) J 



(A.6) 



/ 



o 



T = 



M| x 



(A.7) 



with 






o 



Wan 



(A.i 
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This diagonalization gives the following general result: 
1) massless (normalized) eigenvector 



VSg & 



V3g& 



2V2g 



w 



M\ x 



2) 1st massive neutral boson 
1 



(A.9) 



(4 + S 2 h)(9 2 w + 3^1/8) - J (51 + %)*(<& + 3gl/8y - g 2 w ?> R b\{^ w + ^ 



7,2) 

s) 



(A.10 



with corresponding (not-yet-normalized) eigenvector 



9w 



(8 2 h ~ S 2 r)(9 2 w + 3fl|/8) + V(K + KWw + 3^/8) 2 - g'tv^U^ + W 



7,2) 

s) 



~9sVZ <fo(P H - SI) - m{8l + P H )/8 + v^ + (%)*{& + WW - 9 2 wS 2 R SU^9 2 w + 3£ 



7,2) 

s) 



a/3/2 gw9s&H 



M 2 Z2 



9 2 W (K - S R ) - 3~gl(6l + S%)/8 + ^6% + S' H )\<? W + 3$/8)* - g^tfUWw + 3$) ' 
3) 2nd massive neutral boson 
1 



;1 



(4 + S 2 H )(g 2 w + Sg 2 s /8) + \(S 2 R + 8 2 H Y{g 2 w + 3$/8)* - g 2 w ^\W w + 3g 8 



7.2 s 



(A.ll 



with corresponding (not-yet-normalized) eigenvector 



9w 



(S 2 h - K)(9 2 w + 3^1/8) - v^ + Smg'tv + 3^I/8) 2 - 9 2 W 5' R 5UW W + 3^ 



7,2) 
S) 



~9sVz [gUtl - Sjd - ^W R + s%)/8 - VK + P H n<? w + WW 2 - g'wK^g'w + 3^ 

a/3/2 g w gs$H 



gW\ - 6%) - 3~gl(S R + 8D/8 - y/(P R + P H )*tf w + 3^/8)* - <&<^(4<& + W) 
A special situation where 5 R = 5 2 H has been also discussed in the main text. 



A. 2 Charged gauge boson mixing 

The mixing of charged gauge bosons W^, W R , W% is given by the VEV of A#(4, 1, 3, 1), 
A ff (4, 1,1,3) (see above), $(1, 2, 2, 1), ^(1, 2, 1, 2), <£ 2 (1, 1, 2, 2), 0,(15, 2, 2, 1), $ H (1, 2, 2, 3): 



($) 



Vi/a/2 

^/a/2 



(A.12) 
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(*1> 



u 



i/V2 
u 2 /V2 



(A.13) 



(*2> 



v[ 
v' 



Wx/y/2 

w 2 /V2J ' 

/l \ 
1 10 
2^ 10 
\0 -3/ 



(A.14) 



(A.15) 



(*h) 



/xi/2 \ 

X2/2 

-xi/2 

\ -x 2 /2/ 



with the corresponding covariant derivative 



D^ = d^ - igw[w; L -\® - ig w ^[w; R -i 

^$2 = <9 M $ 2 - »^[W^^]$ 2 - igw*2\W*H^], 



D^ d dp<j> 8 - ig w [W° L - 



-igw< 



m R - h 



~2 



(A.16) 



(A.17) 

(A.18) 
(A.19) 
(A.20) 



± 

L 1 



^0223 = ^0223 ~ WwWl^\ ® 0223 ~^0223 ® [W^y] - l0 W [H^y , 022s]- (A.21) 

From Tr (D^D^cfA follows the squared mass matrix of charged gauge bosons W 

\ 



wl W±. 



9w 



^V 2 + u 2 + X 2 V 2 2 + X \ 2 

V? 2 + Xli 5 2 R /2 + V 2 + w 2 + X 2 



«12 

wj 2 



\ 



u 



12 



™12 



(A.22) 



S 2 H /2 + u 2 + w 2 + Ax 2 ) 



where we have defined 

V 






2 M i + M i 



u 



?/-' 



2 W\ + W 2 2 



(A.23) 



43 



X 



2 xl + xl 



2 ^ XiXj , 



K- 



f,;f,- */*?/?. ..2 



^l/j 



- ',.y 



+ 



UiUj 2 

u i,j = ~ 2~ ' w *,j 



WiWj 



(A.24) 
(A.25) 



In the limit t> 2 , t>' m , w 2 ,, xf <C wf , 5% 6 2 H (as it is considered in the main text), we obtain 
in the leading order approximation 

• 1st massive charged boson's squared mass 






g 2 w(v 2 + u 2 + x 2 + oC ¥ ) 



with the corresponding (normalized) eigenvector 



1 + 0(-„ 



_ 

■52 






°('« 



_ 

'52 



(« 



2nd massive charged boson's squared mass 



M, 



Wi — 9w 



61 + 6 



R~ U H 



+ -J(6l-5] i Y + lQw\ 2 + 0{i 



4 4 

with the corresponding (normalized) eigenvector 



(A.26) 



(A.27) 



(A.28) 



0(v 2 ) ; 2y/2wf 2 + 0( 



2 ^-^+y/(4-^) 2 +16< 

7 ^ ' v^ 



+ C(t; 2 



4 

12 



16< + (51 - 5|,)2 - (51 - 5l)v/(5l-5l)2 + 16^ 
3rd massive charged boson's squared mass 



M k = 9 * w < S « + S » ' 



with the corresponding (normalized) eigenvector 



G(v 2 ) ; 


2V2w 2 12 + 0(v 2 ) ; 


S%-S%-^(5%-6J i r+16wi 2 i c , 2 . 



16< + (51 - 51)2 + (51 - 51)^(51-51)2 + 16™ 



(A.29) 



(A.30) 



(A.31) 



Having these eigenstates, one can straightforwardly construct the orthogonal matrix O 
that diagonalizes the squared mass matrix of charged bosons. The same matrix also 
relates the gauge eigenstates W^, W^, W^ to the mass eigenstates W^, WfS W 2 






/w±\ 



o 



T 






(A.32) 
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In a special case where 5 2 R = 5 2 H = 5 2 , and in the leading order the above expressions 
become quite simple: 



1st massive charged boson's squared mass 



M W =9w[V 2 + u 2 + X 2 + 0(-) 



J 2 ' 



with the corresponding (normalized) eigenvector 



v 2 v 2 v 2 



(A.33) 



(A.34) 



2nd massive charged boson's squared mass 

/5 2 
Ml - ~ 2 ' ' 



l w 1 — 9w I 9 



?L' 



12 



0(v 2 



with the corresponding (normalized) eigenvector 



v 2 I t> 2 I v 2 



(A.35) 



(A.36) 



3nd massive charged boson's squared mass 



M 2 W2 = g 2 w (^ - w\ 2 + 0(v 2 )") 



with the corresponding (normalized) eigenvector 



v 2 1 t> 2 1 v 2 



(A.37) 



(A.38) 



From these eigenstates, one can construct the rotation matrix, 



(wf\ 



O 1 



fwi\ 



which has the following form: 



O 1 



I i o(£) o{$) \ 

0(£) 1/V2 l/v^ 



(A.39) 



(A.40) 
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